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A hydrodynamical model for the FERMI-LAT 7-ray light 
curve of Blazar PKS 1510-089 
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ABSTRACT 

A physical description of the formation and propagation of the working surface inside 
the relativistic jet of the Blazar PKS 1510-089 is used to model its 7-ray variability 
light curve using FERMI-LAT data from 2008 to 2012. The physical model is based 
on conservation laws of mass and momentum at the working surface as explained 
by iMendoza et all (|2009D . The hydrodynamical description of the working surface is 
parametrised by the initial velocity and mass injection rate at the base of the jet. We 
show that periodic variations on the injected velocity profiles are able to account for 
the observed luminosity. With this, we are able to obtain mass ejection rates of the 
central engine which are injected at the base of the jet, and oscillation frequencies of 
the flow, amongst other physical parameters. 
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1 INTRODUCTION 

The study of the 7-ray emission from relativistic jets in AGN 
constitute a modern test for high energy jet-physics. Among 
all types of AGN, Blazars (Blazar class is defined as ra- 
dio loud sources conformed by the BL Lac ob jects and the 
Flat Spectrum Radio Quasa rs -FSRQ, see e.g. Fossati et al.l 
' Il997l : IGhisellini et al.l 1 19981 . and references therein) repre- 
sent the most extreme c lass. They are kno wn to have the 
' most powerful jets (e.g. ILister et al.l [2009) and also show 
a highly variable Spectral Energy Distribu tion (SEP) from 
the radio to the 7-rays wavelengths (see lAbdo et al.ll2010l : 
ID'Ammando et"al]|201ll . and references therein). The emis- 
sion produced by their jets is highly amplified due to Doppler 
boosting effects. This is because the jet of these sources is 
aligned very close to the observer's line of sight. 

Among the Blazar population, the FSRQ PKS 1510-089 
is known to be one of the most powerful objects. Its his- 
torical behaviour shows that it is a highly variable object, 
reaching varia tions as large as ^ 6 mag in the optical B- 
band. (see e.g. lLiller fc L illcr 1975). In the optical, a quasar 
monitoring study done by iLu (|1972T I presents the first op- 
tical light curve of this object. Later, PKS 1510-089 was 
classifi ed as a high polarised quasar by Mo ore fc StockmanI 
|l98ll) . One of the first radio observations done to this ob- 
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ject, at sub-arcseconds r esolution, revealed a clear structure 
of knots in the jet fsee lO'Dea et al]|l988h . These authors 
also find that the observed hot spots had velocities ~ 0.93c. 
Recently, it has been found that PKS 1510-089 has a highly 
collimated relativistic jet with apparent superluminal speeds 
going from 20 c to 46c and with a s emi-angle aperture for the 
jet of ~ 0.2° (|jorstad et al.ll2~005l ). Since the angle between 
the relativistic jet and the observer's line of sight is 1.4° - 
3°, the jet almost coincides with the obse rver's line of sight 
i|Homan et al.ll2002i ; iMarscher et al.ll2010l ). 

PKS 1510-089 was one of th e 7-ra y sources de- 
tected by EGRET (jHartman et all Il999h . Recently, it 
began t o be monitored a t high energies with A GILE 
(see e.g. iPucella et all 120081 : ID'Ammando et all 120081) and 
later also by FERMI-LAT and AGILE ilTramacereH2o"o"8l ; 
ICiprini fc Corbell 120091 : ID'Ammando et all 120091 ) It has 
also been studied with MAGIC and HESS l|Cortinal 
|2012| ; I Wagner etahl 120101 ). The most prominent out- 
bursts displayed by PKS 1 510-089 were observed in 2008 
l|Kataoka et al.ll200Sft . 2 009 (ICiprini fc Corbelll2009l ) and the 
most recent in 2011 (see lOrienti et al.ll2012r )~ Therefore, the 
high activity observed in this source, turns it into an ideal 
target for the physical study of its highly relativistic jet. 

Precise models for the light curve produced by the out- 
burst and flares from Blazars are not done using directly 
the data variations observed in different wavelengths. In- 
stead, models are applied to explain the behaviour of the 
Spectral Energy Distribu tion (SED) (e.g. lAbdo et al1l2010l ; 
ID'Ammando et "aT]|201ll ). Direct understanding of the light 
curve requires a precise knowledge of the hydrodynamical 
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behaviour of the relativistic flow. iMendoza ei~ai1 (|2009h have 
constructed a hydrodynamical model of the motion of a 
working surface inside a relativistic jet which is able to fit the 
observed light curves of long Gamma- Ray Bursts (lGRB's). 
Since the jets in Blazars are highly relativistic and their jet 
is nearly pointing towards the observer, similar to the jets 
observed in lGRB's, the physical ingredients of both phe- 
nomena can be considered the same but occurring at differ- 
ent phy sical scales of energy, sizes, masses, accretion rates, 
etc. (cf. lMirabel fc Rodriguezll2002h . 

The Blazar PKS 1510-089 is of tremendous importance 
since it exhibits extreme relativistic motions. As such, its 
energy curve must present luminosity variations and peri- 
ods of extreme activity displayed as outbursts that, when 
physically modelled, can yield a better understanding of the 
physical parameters associated to the mechanism producing 
the observed luminosity. 

In this letter, we assume that the mechanism producing 
the observed light curve in the 1GRB is exactly the same that 
produces the variable light curve of the Blazar PKS 1510- 
08 9. We thus apply the h ydrodynamical jet model presented 
in IMendoza et all |2009) to the light curve variations dis- 
played by the Blazar PKS 1510-089 in the 7-ray domain, 
through observations done with the FERMI-LAT telescope. 

The letter is organised as follows. In Section [2] we 
explain in general terms the data reduction process. In 
Section [3] we describe the characteristics of our hydrody- 
namic model. The fit done to the data with the hydro- 
dynamic model is explained in Section [4] The results of 
our fits and the discussion of the main physical parame- 
ters obtained in the modelling are presented in Section [5] 
Throughout this paper we use a standard cosmology with 
Hp = 71kms~ 1 Mp c" 1 , Q m = 0.27 and fi A = 0.73 (see e.g 
iKataoka et al.ll2008l . and references therein). 



2 FERMI-LAT DATA 

The gamma-ray fluxes were obtained in the range of 0.2 to 
300 GeV using the public database of FERMI-LAT from 
2008 August 08 to 2012 May 28. The data were reduced 
with the FERMI science tool package (see e.g. lAtwood et all 
2009) in the same energy range, taking into account the dif- 
fuse galactic background radiation, the instrument response 
matrix p7v6, and considering a zenith angle < 105°. We also 
calculated the active time of the detector and the PSF. The 
7-ray light curve (LC) was constructed modelling the flux 
with a power law of the form dN/dE = N (E/E )~<. The 
fluxes and errors obtained with this package are given in 
photons x cm _2 s _1 . For further physical interpretation of 
the data, we have converted these fluxes and errors to MeV 
cm~ 2 s" 1 . 

The photons considered for analysis were taken from a 
region centred on the coordinates of PKS 1510-089 with a 
radius of 15°. Figure[T]shows the 7-ray LC, with a bin size of 
1 day. We chose these bins, since the errors are larger using 
shorter bin sizes, complicating the analysis of the data and 
because particular outbursts can be adequately resolved. 

From Figure [1] it follows that the source displayed the 
historical maximum outburst in MJD 5585 1, correspond- 
ing to 2011 October 17 and reported by iHungwe et al.l 
|201lT ). Another important outburst occurred in MJD 



54899 (2009 March 9) and was observed with AGILE 
l|D'Ammando et alj [2009 ) . Several flares or outbursts can 
be observed in the LC. The most relevant events occurred 
in MJD 54717 (2008 September 8), M JD 54843 (2009 Jan - 
uary 12), MJD 55200 (2010 January 4. lBemtez et alj|201ll ). 
MJD 55730 (2011 June 18), and MJD 55954 (2012 Jan- 
uary 28). This last e vent was also observed by A GILE 
JVerrecchia et alj|2012|) an d MAGIC (|Cortinal \20vh . Note 
that iMarscher et al.l (|2010T l reports extra flares < 200MeV 
during the period 54850 - 54950 MJD, which are not seen 
in our > 200MeV selection. 



3 A HYDRODYNAMICAL MODEL FOR THE 
LIGHT CURVE OF PKS 1510-089 

The formation of internal shock waves on a relativistic jet are 
commonly explained by different mechanisms, such as the 
interaction of the jet with inhomogeneities of the surround- 
ing medium, the bending of jets and time fluctuations in the 
parameters of the ejection (see e.g. iRees fc Meszaros! Il99~ " 



Mcndoza & Longair 2002; Jamil et al. 2008: Mcndoza et al 



2009). In particular, the model by IMendoza et al.1 |2009l ') is 



a hydrodynamical description that can be applied to shock 
waves inside relativistic jets. This semi-analytical model de- 
scribes the formation of a working surface inside a hydrody- 
namical jet due to periodic variations of the injected flow. 
When fast flow overtakes slow flow, an initial discontinuity 
is formed and a working surface (two shock waves separated 
by a contact discontinuity) is produced. The working sur- 
face travels along the jet and radiates away kinetic energy. 
IMendoza et al.l (|2009T ) assumed that the efficiency convert- 
ing factor is ~ 1 and that it is mostly emitted in the 7-ray 
band. As explained in Section [TJ the Blazar PKS 1510-089 
behaves as a n scaled typical 1GRB and as such, the hypoth- 
esis used by IMendoza et al.1 (2009) can be extended to this 
particular object. As we will discuss in section O this as- 
sumption is coherent with the physical prop erties found from 
the model. Following IMendoza etafl (2009), we assume that 
flow is injected at the base of the jet with a periodic velocity 
given by 



■j(t) = v + erf 



(1) 



where r is the time in the rest frame of the source, the ve- 
locity vo is the "average" velocity of the flow inside the jet, 
and cj is the oscillation frequency. The positive constant pa- 
rameter rf is chosen in such a way that oscillations of the 
flow are small so that the bulk velocity v(t) of the flow does 
not exceed the velocity of light c. The mass ejection rate 
m(r) from the central engine which is injected at the base 
of the jet is assumed constant through a particular outburst 
event, but is allowed to vary from one outburst to another. 
The radiated energy of the flow as a function of time is 
calculated as the difference between the total energy Eo in- 
jected at the base of the jet and the kinetic energy inside 
the working surface E ws - The luminosity L is thus calcu- 
lated as_th^_derivatwe ofthis radiated energy with respect to 
time. IMendoza et al.l (|2009T ) showed two ways of calculating 
this luminosity curve. The first method consisted in a semi- 
analytical procedure and the second is performed with a full 
hydrodynamical numerical model. The authors showed that 
the semi-analytical model is in good agreement with the full 
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Figure 1. Fermi-LAT light curve of Blazar PKS 1510-089 (from 0.2 to 300 GeV) obtained from 2008 August to 2012 May. The outburst 
identification number (ID#) labelled in the figure stands for the different flares selected in our work (see text). The 3c noise level is 
represented by the red horizontal dashed line. 



numerical simulation, and as such we model the light curve 
of PKS 1510-089 using their semi-analytical approach. 

The semi-analytical approach is based on the assump- 
tion that equation fl]) is valid and as such, one needs to 
know (or find through fits to observational data) the val- 
ues of vo, rj 2 , u) and m. Furthermore, in the model, the 
mass ejection rate rh enters in the description of the prob- 
lem since the luminosity L oc the 2 . The value of vo can be 
safely taken as the average velocity of the flow of the jet, 
which must come from observat ional data (for this particular 
source [D'Ammando et al]|2008l . reports a value r(i>o) = 18). 
Since the value of rj 2 has to be small, one starts with a value 
of rj 1 such that the bulk velocity of the flow v ~ 0.1 X vq 
and one allows it to vary in such a way that the Lorentz 
factor r(«) of the full velocity v never exceeds the upper 
limit of ~ 200, according to the standard expectations of an 
very large upper limit for the Lore ntz factor associated with 
Blazars (see e.g. iLister et al.| [2009). With this, the model is 
left with three parameters (rj 2 , rh, u) that fit the observed 
features in the light curve. 



4 MODELLING THE 7-RAY LIGHT CURVE 

To model the LC of Figure [TJ we have selected the most 
conspicuous flares and outbursts. The criterion used consists 
of selecting only those flares or outbursts that are beyond 
3cr noise level according to the errors shown in the LC. By 
doing so, it turns out that 38 relevant peaks were chosen for 
our fitting. 

As explained in section the model has four free pa- 
rameters. The velocity parameter Vq is fitted by assuming 
an average bulk velocity for the jet such that its Lorentz 
factor is T(vp) = 18 a ccording to the observations of 
iD'Ammando et all |201ll ). For a given outburst, we obtain 
the parameters rh and w by doing a linear regression analy- 
sis on_tfte_d^ta L _u^ng_the model of Section O since according 
to iMendoza et all (|2009r ). L oc mc and t oc w . The ex- 
act value of rj 2 is obtained by trial and error until the best 
fit (to within 10% of accuracy in uj and m) is obtained. To 
calculate the Flux F we divide the obtained Luminosity L 
by AitD 2 L , where Dl is the luminosity distance, which for 



this particular case is Dl = 1919 Mpc. The results of these 
fits are shown in Figure [2] The obtained values of the phys- 
ical parameters of the model for any particular outburst are 
presented in Table [T] 

There are a certain subclass of outbursts that we do 
not model. These outbursts, labelled 8, 10, 20, 27 and 32 in 
Figure [1] do not have enough data to allow us an accurate 
modelling. The outburst labelled 11 seems to have a fall that 
develops into a constant value before reaching an expected 
minimum and no data points further, so it seems incomplete. 
The outbursts 34 and 35 have large errors which also makes 
the modelling not accurate. 



5 DISCUSSION 

Table [T] shows the results of the fitting performed on 29 
outbursts of the LC of PKS 1510-089. The obtained values 
of the parameter rf are chosen in such a way that the errors 
in lu and rh do not exceed 10%. We have also included the 
maximum and minimum Lorentz factors, P m ax and r m i n , 
obtained for the bulk velocity of the flow. 

We have modelled the Light Curve of Blazar PKS 1510- 
089 for almost 4 year s using the hydrodynamical model of 
IMendoza et al.l (|2009l ). The modelling was performed by as- 
suming a periodic velocity injection mechanism at the base 
of the relativistic jet that leads to the formation of a working 
surface and is capable of loosing energy as it t ravels along the 
jet. As explained in section [3j the model by IMendoza et all 
was constructed to deal with light curves of 1GRB. However, 
the Blazar PKS 1510-089 has many physical characteristics 
to be considered a geometrical large scaled version of a 1GRB 
since it has a highly relativistic jet that points towards the 
observer. In fact, the results presented in Section[S]and Ta- 
ble [1] are in good agreement with the expectations of the 
physical values that a Blazar jet can have, except for a few 
large (> 150) bulk Lorentz factors founded for peaks 13, 
16, 25, 27 and 30. These somewhat large Lorentz factors in 
the velocity oscillation show another close similarity with 
lGRB's. 

The resolution of one day in the FERMI-LAT data 
allowed us to study in great detail the outbursts of PKS 
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Figure 2. In each panel, the continuous red curve re presents the fitting don e to the light curve of PKS 1580-089 with the semi-analytical 
model of internal shock waves (working surfaces) bv lMendoza et al ] 42009T ). Blue horizontal dotted-lines in all panels show the 1 a and 
3 a noise levels. Top-left panel shows variations from 2008 September to 2009 February. Top-right shows variations from 2009 March 
to 2009 May, where a couple of large outburst occurred. At the end of these events, a more quiescent state was observed, although a 
series of successive small flares can be identified. Central left-panel shows a zoom of the peaks 12 and 13. Central right-panel shows 
the variations observed from 2009 December to 2010 March. Bottom left- panel shows recent variations occurred from 2011 July to 2012 
March. Finally, bottom-right panel shows a zoom of the October 2011 outburst. This outburst is ~ three times more luminous than the 
previously one observed in 2009 March. Up to now, this is the most violent outburst observed in the 7-ray waveband by FERMI. 



1510-089. Some of the outbursts reported in this letter, cor- 
responding to early 2012, were als o reported by AGILE 
l|Lucarelli et alj|2012h and MAGIC l|Cortinj|20l"3 ). 

The range of parameters as presented in Table [T] i.e. 
rh ~ 10 -7 - 10 _10 M©yr _1 , ui' 1 ~ 20 - 70min and varia- 
tions of the Lorentz factor T between 20—18 0, denote a scal- 
ing be tween the 1GRB counterparts found in lMendoza et all 
(2009) for which rh ~ 10" 1 - 10" 2 M©s _1 , u' 1 ~ 10s and 
r ~ 50 — 500. The inferred high relativistic Lorentz fac- 
tors associated to the motion of the bulk velocity of the 
flow inside the jet of PKS 1510-089 makes it an ideal can- 
didate for the a pplic ation of the hydrodynamical model of 
iMendoza et alj (|2009l ) . This is why that physical model can 



be applied naturally to 1GRB and in this particular case 
to the extreme relativistic motion of the jet in the Blazar 
PKS 1510-089. The variations of the injected flow at the 
base of the jet cause the formation of working surfaces that 
produce flares or bursts of 7-rays in the structure of the jet. 
The physical mechanism producing the oscillations of the 
input flow, which allows fast fluid to overtake the slow one, 
leading to the formation of working surfaces, is beyond the 
scope of this letter. However, steady flow deviations and os- 
cillations in such complicated phenomena are expected since 
the accretion-ejection mechanism associated to a particular 
object is not necessarily expected to be of constant velocity 
and mass accretion-ejection rates. 
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Table 1. Different physical quantities obtained for the outbursts 
modelled in this work. The mean Lorentz factor of the bulk ve- 
locity of the flow was assumed to be 18. The first three columns 
from left to right are the period, numeric identification of the 
outburst (ID #) and the date corresponding to the maximum lu- 
minosity for a particular outburst. Columns four and seven are 
the obtained values for the parameters r/ 2 (measured in units of 
the speed of light c) and the inverse frequency u)~ 1 relevant to 
the particular variational model of equation JTJ. The errors asso- 
ciated with these parameters have an accuracy of 10%. Column 
eight represents the mass injection rate fa of the flow at the base 
of the jet. Columns five and six correspond to the minimum and 
maximum values of the Lorentz factor of the bulk flow for each 
particular outburst. 



It is important to note that the assumption of seeing 
a Blazar as a scaled version of a 1GRB is not new. In an 
early attempt for finding a unified model o f jet and central- 
engine power, iMirabel fc Rodriguez! (|2002l ) made this iden- 
tification. The more relativistic a Blazar jet is, the more it 
will resemble a 1GRB. The idea of having a unified physical 
model for all types of astrophysical jets was first suggested 
by the pionee ring works for the astro physical sca ling laws of 
black holes bv lSams et ahl (|l996l ) and lReesI ljl998h . The work 
presented in this letter further strengths arguments about a 
unified picture of all astrophysical relativistic jets. 

PKS 1510-089 result e d to b e an ideal target to test the 
model bv lMendoza et al.l (2009) since it closely resembles a 
1GRB in some of its outbursts. Future tests of the model 
have to be done with a wide variety of Light Curves from a 
large collection of Blazars and micro-quasars. 
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